Diagonal-unitary 2-designs and their implementations by quantum circuits 
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We study random diagonal-unitary matrices, an ensemble of unitary matrices diagonal in a given 
basis with randomly distributed phases for their eigenvalues. We investigate how to efficiently 
implement random diagonal-unitary matrices by quantum circuits and introduce diagonal-unitary 
t-designs that simulate diagonal-unitary matrices up to the t-th order moments. We present two 
efficient implementations of diagonal-unitary 2-designs in the computational basis by using quantum 
circuits. One is composed of single-qubit diagonal gates and controlled-phase gates with random- 
ized phases, which achieves an exact diagonal-unitary 2-design after applying 0(N 2 ) gates. If the 
controlled-Z gates are used instead of controlled-phase gates, the circuit achieves an approximate 
2-design and it is necessary to apply 0(N 2 (N + logl/e)) gates. We also provide applications of 
the diagonal-unitary 2-designs implemented by the quantum circuits, a generation of a 2-design of 
random states and a quantum simulator of canonical states of classical Hamiltonians. 

PACS numbers: 03.67.Ac, 03.67.Bg, 05.67.-a 



I. INTRODUCTION 

Random unitary matrices are an ensemble of unitary 
matrices uniformly distributed in terms of the Haar mea- 
sure on the unitary group. They have been studied in 
the context of random matrix theories [l[ and arc re- 
cently used for applications in quantum information pro- 
cessing [2I-Q. Since random unitary matrices provide a 
method of random sampling of quantum states, they are 
also used for investigating typical properties of phy sical 
systems when derails of the system are unknown [l3,|TH. 
An ensemble of quantum states generated by random uni- 
tary matrices are known as random states and have been 
particularly investigated for studying typical properties 
of quantum states [12T[26j . 

Efficient implementations of random unitary matrices 
are necessary to exploit them for quantum information 
processing [27l . [28j . In particular, the implementations 
by quantum circuits are intensely studied [3, H, I29U3H . 
It has been shown that the ensemble of unitary matri- 
ces simulating up to the t-th order of statistical moments 
of random unitary matrices, which is referred to as ran- 
dom unitary t- design, can be approximately achieved by 
quantum circuits. The number of elementary gates in the 
circuits scales with 0(Nt 4 logt(Nt+\og 1/e)), where N is 
the system size and e is an allowed error for the approx- 
imation |3lj |. This result implies that an approximate 
t-design of random states can be efficiently generated by 
the quantum circuit. 

Besides the applications for quantum information pro- 
cessing, the random unitary matrices provide a method 
to investigate typical dynamics in uniform systems. How- 
ever, it is often the case that we would like to investigate 
dynamics under some constraints of systems. This leads 
to an idea to investigate an ensemble of unitary matri- 
ces with some constraints. For instance, when the sys- 
tem has a symmetry, the Hilbcrt space is decomposed 
into irreducible subspaces and it is natural to consider 



dynamics in each subspace separately. Here, we con- 
sider another constraint of dynamics originated from the 
equation of motion, namely, the Schrddinger equation. 
When the Hamiltonian is time-independent, the dynam- 
ics is represented by a diagonal unitary matrix in the 
eigenbasis of a given Hamiltonian. Thus, it is worth in- 
vestigating random unitary mateices diagonal in a fixed 
basis, which we refer to as random diagonal-unitary ma- 
trices. When random diagonal-unitary matrices are ap- 
plied to a fixed state, we obtain phase-random states, an 
ensemble of states with randomized phases [32j . Phase- 
random states have a close connection to the study of 
time-average properties in quantum statistical mechan- 
ics. 

In this paper, similarly to random unitary matrices, 
we consider implementations of random diagonal-unitary 
matrices by quantum circuits and their applications. To 
this end, we introduce diagonal-unitary t-designs. We 
show that a diagonal-unitary 2-design is exactly achiev- 
able by a phase-random circuit composed of 0(N 2 ) di- 
agonal gates if we are allowed to use controlled-phase 
gates with random phases in addition to single qubit 
phase gates. On the other hand, if we use controlled-Z 
gates instead of controllcd-phasc gates, an approximate 
2-design is achieved after applying 0(N 2 (N + logl/e)) 
diagonal gates. These results show that random param- 
eters in genuine two-qubit gates enhance the ability of 
randomizing phases. We also provide two applications 
of diagonal-unitary 2-dcsigns implemented by the phase- 
random circuit. First, we show that, by combining the 
phase-random circuit with a simple classical procedure 
requiring 0(N) random bits, we can generate an exact 
2-design of random states. Our method of generating 
the 2-design of random states is simpler than previously 
known implementations, particularly from an experimen- 
tal point of view. Second, we propose a quantum simula- 
tor of canonical states of classical Hamiltonians by using 
the phase-random circuit where we do not have to pre- 
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pare the Hamiltonian itself. 

This paper is organized as follows. In Sec. [H] we give 
definitions of random unitary matrices, random diagonal- 
unitary matrices, and corresponding designs, a unitary 
t-design and a diagonal-unitary t-design. In Sec. IIII1 we 
present our two main theorems on the convergence of 
the phase-random circuits and provide the applications. 
We present proofs of the first theorem and the second 
theorem in Sec .11 VI and in Sec.[V] respectively. Finally, we 
summarize and present concluding remarks in Scction lVTl 



II. RANDOM MATRICES AND DESIGNS 

We present basic definitions of random (diagonal- 
unitary matrices, (phase)-random states and the corre- 
sponding designs. 



A. Random unitary matrices and random 
diagonal-unitary matrices 

Definition 1 (Random unitary matrices) 

Random unitary matrices Uuaar are the ensemble 
of unitary matrices uniformly distributed in terms of the 
Haar measure. 

Definition 2 (Random states) Random states 
Yrandom are the ensemble of states {U^ |^ , )}c/ M eWHaar f° r 
any fixed state \^>) £ H. 

Note that the distribution of random states is indepen- 
dent of the choice of the state |^) since the Haar measure 
is the unitarily invariant. 

We also define random diagonal-unitary matrices, 
which are understood as random unitary matrices with 
a fixed basis. 

Definition 3 (Random diagonal-unitary matrices) 

Random diagonal-unitary matrices in an orthonormal 
basis {|it n )}, Wdiag({|w„)}), are an ensemble of diagonal 

unitary matrices of the form U v = Yln=i e%iPn \ u n)( u n\ 
where the phases ip n are uniformly distributed ac- 
cording to the normalized Lebesgue measure dip = 
dpi ■■■dp d /(2n) d on [0, 27r) d . 

Definition 4 (Phase-random states) For a given 
state = J2 n T nZ %u)n \ u n) where r n > 0,uj n € [0, 27r) for 
all n, phase-random states Tph a se({?"n, |w n ) n }) are an en- 
semble of states {U v \^)}u v eU dia . B ({\u„)}), which is equiv- 
alent to {J2 n e lVn r n |Un)} v „e[0,27r)- 

Random diagonal-unitary matrices are closely related 
to the study of the typical properties of time evolutions 
by time-independent Hamiltonians. Under the assump- 
tion that the phases of the expansion coefficients in the 
cigenbasis of a Hamiltonian are randomly distributed af- 
ter sufficiently long time, investigations of the statisti- 
cal properties of random diagonal-unitary matrices are 



equivalent to the studies of typical properties of such 
Hamiltonian dynamics. Phase-random states are intro- 
duced for this aim and it has been shown that canonical 
distributions are realized in subsystems when an initial 
state and the eigenstates of a Hamiltonian satisfy a trade- 
off relation ^2l K52| . 



B. Designs 

In general, a design of another ensemble is an ensemble 
simulating average properties of the original ensemble. A 
unitary t-design is an ensemble of unitary matrices that 
simulates up to the t-th order of moments of random 
unitary matrices on average. 

Definition 5 (Unitary t-designs) Let Wnaar be ran- 
dom unitary matrices. A unitary t-design U^L^ is a dis- 
tribution of unitary matrices ti^ = {U^d^ such that 



u? ® (utr'dn, 



where W is an Hcrmitc conjugate of U. 

Corresponding designs for random states are re- 
ferred to spherical t-designs [3314351 ] . complex-protective 
t-designs (3(1 l37j or quantum state t-designs [2!|, which 
are defined as follows. 

Definition 6 (Complex projective t-designs) 

Let T ranc iom be random states. A complex projective 
t-design Tcp is a distribution of quantum states {\ipfi)}dfi 
satisfying 



|lM6Tra„do 



When the distribution of the designs is discrete, the 
integral is replaced by a summation over a probability 
distribution, which is sometimes referred to as weighted t- 
designs [38j . By definition, a complex projective t-design 
cannot be distinguished from random states if we have 
only t copies of states. 

In analogy with unitary t-designs and complex pro- 
jective t-designs, we define designs for random diagonal- 
unitary matrices and phase-random states. 

Definition 7 (Diagonal-unitary t-designs) Let 

^diag({|w ra )}) be random diagonal-unitary matrices in a 
basis {|it„)}. A diagonal-unitary t-design ^[ ag {{\un)}) 
is a distribution of unitary matrices {t/^}^ diagonal in 
the basis {|u„)}, which satisfies that 



us* ® (uir'dip 



uf ® (uir'dp,. 
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Definition 8 (Toric t-designs) Let 

Yphasedjm be phase-random states. A toric 

t-design T^* r ({r„, |u n )}) is a distribution of states 
{\ipn)}dn satisfying that 

/ , ,\<P>p){^ v \ 9t dtp = / {^(ip^dcp. 

Since the parameter space of phase-random states is a 
c?-dimensional torus, we call this type of designs as toric 
t-designs. 

The state designs can be generated by applying unitary 
designs to a state. Inversely, if a set of unitary matrices 
(diagonal-unitary matrices in the basis generates 
a complex projective t-design (a toric t-design) for any 
state \^f), then the set is a unitary t-design (a diagonal- 
unitary t-design). 

In this paper, for simplicity, we denote by E any ex- 
pectations over a probability distribution. When it is 
necessary, we explicitly write the space taken over for 
the expectation, e.g., E^ Haar an d ^x< 4) ■ Note that the ex- 
pectation of unitary matrices is not necessarily unitary 
and the expectation of states is generally a mixed state. 



C. An example of diagonal-unitary t-designs 

We show a simple example of a diagonal-unitary t- 
design where the phases are randomly chosen from a 
discrete set. This example shows that for achieving 
diagonal-unitary t-designs, continuous random parame- 
ters such as a e [0, 2ir) are not necessary and we can use 
(t + l)-valucd discrete random parameters instead. 

Proposition 1 A set of unitary matrices fit = 
{J2n eZrj,n \ u n) ( u n\}{ui n }, where (f> n is randomly chosen 



2nk 
t+1 

basis {!««)}• 



from {|rf }fe=o,i,--- .*■ is a diagonal-unitary t-design in the 



Proof 1 For U+ = £„ e*" \u n )(u n \, Uf l ® 17+®* is cal- 
culated to 



U*t®Ul»* = £ exp^(^ fc -0 m J] 

ni 1 mi ,mt — 1 k— 1 

I ^ni ' ' ' V>rit %i ' ' ' 'Unit ) (^^l ' ' ' ^Tlt ^ l ' ' ' ^ | 

For and fl t , the expectation of an operator X is 

taken over G [0, 2ir) and <\)\ <E {fr|}fc=o,i,— ,t for all 
* = !,••■ , d, respectively, namely, 



E,, t) [X] 



1 



Xd(f>i ■ ■ ■ d(j>d, 



En t [X] = ( TTT ) d 



E - E 



X. 



The equation E Qt [U®* <g> f/t®*] = E (t) [C/®< ® C/t®*] fol_ 

diag 

lows from an identity that 



(2 



/ exp[i^(<^> 

n k 4* in k 

)}d<f>i ■ ■ ■ dcfid 

n ' ^° fe=i 



(j^r J2 E exp[i^(^-^J]. 



Generation of the matrices of diagonal-unitary t- 
designs based on Proposition [1] requires global random- 
izations over the Hilbert space of n qubits. For implene- 
tating these matrices by using quantum circuits, we need 
to decompose diagonal-unitary t-designs into local uni- 
tary operations, namely, one- and two-qubit gates, in an 
efficient way. This is the main concern of this paper. 



D. e-approximate t-designs 

An e-approximatc t-design is an ensemble that approx- 
imates the t-design within an error e [H . For unitary de- 
signs and diagonal-unitary designs, we use the diamond 
norm to evaluate the difference between two matrices. 
For a superoperator £ on H, the diamond norm is de- 
fined by 



sup sup 

d X^O 



\(S ®id d )X\i 



where id<j is an identity operator on another d- 
dimcnsional Hilbert space H! and X is any positive oper- 
ator on % <8> H! [3!| . The diamond norm gives the prob- 
ability of distinguishing two operations if we are allowed 
to use an auxiliary system with any dimension. 

Definition 9 Let U^> be a unitary t-design or a 
diagonal-unitary t-design. An e-approximatc t-designs 
of WW is a set of unitary operators U^' e ' that satisfies 



< e. 



For state designs such as complex projective and toric 
t-designs, we define an e-approximatc t-designs in terms 
of the trace norm. 

Definition 10 Let T^) be a complex projective t- 
design or a toric t-design. An e-approximate t-designs 
of T^*) is a set of states T^'^ that satisfies 



||E T(t) [|*}(*| w ]-E T (^[|*}(*l 
where \\X\\ = Ti\X\ is the trace norm. 



i < c, 



Since the t-designs are sufficient to perform most of 
known quantum information tasks using random uni- 
tary matrices or random states [40| . they have been 
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FIG. 1: A phase-random circuit. The vertical line denotes 
a two-qubit gate W^ 1 *'""^ randomly selected from a diagonal 
two-qubit gate set Wdi ag on randomly chosen a pair of it-th 
and jt-th qubits. 



intensely studied. Many exact or approximate com- 
plex projective i-designs are known for iV-qubit systems 
HI) HE SH 51] • The construction of unitary t-designs 
by using quantum circuits has been also studied [H, [H, [29i — 

Little is known about how to construct random 
diagonal-unitary ^-designs and toric i-designs by quan- 
tum circuits. In Ref. [321 ] . it is shown that a quantum 
circuit composed of controlled- Z and single qubit phase 
gates generates an ensemble of states of which average 
amount of entanglement is same as that of phase-random 
states measured by a linear entropy of a reduced density 
matrix. However, it had not known whether the circuit 
achieves a diagonal- unitary i-design or not. 



III. MAIN RESULTS AND APPLICATIONS 

In this section, we introduce phase-random circuits 
proposed in [32j and state our main results, i.e., the 
phase-random circuits can achieve diagonal-unitary 2- 
designs. Then, we provide two applications of the phase- 
random circuits. 

Note that diagonal-unitary t-designs for t > 2 can 
generate large amount of entanglement [32|, which im- 
plies that diagonal-unitary t-designs for t > 2 cannot 
be achieved if we use only one-qubit gates since one- 
qubit gates cannot increase entanglement. We also note 
that diagonal-unitary i-designs for general t cannot be 
achieved by applying random diagonal one- and two- 
qubit gates since the commutability of gates leads to the 
lack of parameters. Therefore, we concentrate in the case 
of t = 1, 



A. Phase-random circuit 

We denote by \n) the computational basis where n 
is a binary representation of n (n = 1, • • • ,2^). We 
investigate achievability of a diagonal-unitary i-design 
in the computational basis W^„({|n}}) by using phase- 
random circuits presented in the previous work of the 
authors (32| . For implementations of a diagonal-unitary 
t-design in a general basis, it is sufficient to apply a uni- 



tary operation transforming the computational basis to 
the basis before and after applying the phase-random cir- 
cuit. 

A phase-random circuit consists of T diagonal two- 
qubit unitary gates shown in Fig.[TJ For the t-ih gate, we 
select two different numbers (it,jt) from {1,2,- •• , N}, 
as well as a two-qubit gate Wt randomly from a given 
set of diagonal two-qubit gates Wdiag- We apply the 
two-qubit gate W t on it-th. and j't-th qubits. An in- 
stance of the circuit is then specified by a set of pa- 
rameters, Ct ■— {it, jt, Wt}J = i, and the unitary op- 
eration represented by the circuit is given by Ut = 



W, 



(iT,jT)yy(iT-ljT-l) 



T-l 



(*i>Ji) 



where W, 



(it,jt) 



acts on 



it-th. and jt-th qubits. Thus a phase-random circuit con- 
sisting of T two-qubit gates is denoted by a set of the 
unitary operations {Ut}c t - 



B. Main results 

The main result of this paper is that, if we choose an 
appropriate two-qubit diagonal gate set Wdiag, the phase- 
random circuit achieves a diagonal-unitary 2-design in 
the polynomial number of the gates. The necessary num- 
ber of the gates depends on the choice of the gate set 

Wdiag- 

We consider two gate sets. First, we study a gate 
set given by w£f g = {diag(l,e ia ) ® diag(l,e ip ) ■ 
CP(l)} a ^ l£{0 ^ m where CP( 7 ) = diag(l,l,l,e*r) 
is a controlled-phase gate with a phase 7. We refer to 
the phase-random circuit consisting of this gate set as 
a CP phase-random circuit, where the set of parameters 
is given by C^ p = {at, Pt, 7t}t=i- Then, we obtain the 
following theorem: 

Theorem 1 The CP phase-random circuit is a 
diagonal-unitary 2-design in the computational basis if 
two-qubit gates randomly chosen from a gate set W^^ 
are applied on all the different pairs of qubits. Thus, the 
number of the required gates is given by T = ^ ■ 

Since all gates in the CP phase-random circuit are com- 
mutable, they can be applied simultaneously in a prac- 
tical implementation so that the depth of the circuit is 
0(1). 

Next, we deal with a gate set given by W^^ = 
{diag(l,e ia ) ® diag(l,e if} ) ■ CZ} a [j£{0 ^ where 
CZ := diag(l, 1, 1, — 1) is a controllcd-Z gate. In 
this case, we have to choose (it,jt) randomly at each 
time, so that the set of parameters is given by Cq! z = 
{it,jt,°tt,/3t}t=i- We call the corresponding phase- 
random circuit as CZ phase-random circuit. Due to the 
fact that a controlled-Z gate has no parameter, the circuit 
converges slowly as stated in Theorem [2j 

Theorem 2 The CZ phase-random circuit {Ut\c° z 
consisting of T two-qubit gates randomly chosen from 



5 



a gate set W^f g applied on a pair of randomly chosen 
qubits is an e-approximate diagonal-unitary 2-design if 
T > T conv (e) where 

N 3 N 2 

— log 2 + — log e- 1 + 0(N 2 ) < T conv {e) 

< 3N 3 \og2 + N 2 log e" 1 +0(N 2 ). 

Therefore, the CZ phase-random circuit is an e- 
approximate diagonal-unitary 2-design after applying 
0(N 2 (N + log e" 1 )) two-qubit gates. 

In a practical sense, the CZ phase-random circuit has 
disadvantagescomparing to the CP phase-random cir- 
cuit since it cannot achieve an exact diagonal-unitary 
2-design. Moreover, unlike the CP case, the gates in 
the circuit cannot be applied simultaneously since the 
dynamics should be stochastic by choosing i and j ran- 
domly for each gate in order to sufficiently randomize a 
state, which will be seen in the proof of Theorem^ How- 
ever, we present Theorem [2] since the difference between 
the CZ and the CP phase-random circuits show that the 
additional parameter of the controlled-phase gates dra- 
matically improve the ability of randomization, which is 
interesting in the theoretical point of view. 

Note that, for both of phase-random circuits, it suf- 
fices to take the phases a, j3 and 7 from a discrete set 
{0, =^-, instead of a continuous set [0, 2ir) as shown in 
the example with t = 2 of the diagonal-unitary t-designs 
presented in Subsection III CI 

Finally, we emphasize that the phase-random circuits 
are expected to be easily implemented in experiments 
since they are composed only of up to two-qubit gates 
diagonal in the computational basis. Furthermore, in 
the case of the CP phase-random circuit, all gates can be 
applied simultaneously, which significantly simplifies the 
experimental implementation. 



C. Applications of phase-random circuits 

We show applications of the phase-random circuits. 
We propose applications for generating a complex pro- 
jective 2-design and for quantum simulators of canonical 
states. 



1. Generating a complex projective 2-design 

The phase-random circuit can generate toric 2-designs 
exactly for the CP case or approximately for the CZ case. 
By combining the phase-random circuits with an extra 
classical procedure, we can also obtain a complex projec- 
tive 2-dcsign. To show this, we first consider the differ- 
ence between complex projective and toric 2-designs. 

The expectation of states for complex projective t- 



designs is calculated to 



®«1 



|Vv)eT randl 

— nw 

j "sym J 



iw>«vr^ 



(i) 



where niy m is a projector onto the symmetric subspace 
in "H®* and d syra is the dimension of the symmetric sub- 
space. Equation (fTJ) is obtained from Schur's lemma [42j . 

On the other hand, the expectation of the state for 
toric i-designs is not a projector IIsy m . For instance, in 
iV-qubit systems, the expectation of the state for the toric 
2-designs T^* r ({r n , \u n )} n ) is calculated to 



E T (2)[|V>) (ip\® 2 } = \ u nU n ){u n u n \ 



2 E r n r ' 



\u„U m ) + \u m U n ) (u„U m \ + (u m U n \ 



V2 



^/2 



(2) 



which is not proportional to the projector II 

To close the gap, we consider a particular toric 2- 
design given by Y^ r ({2 -Ar / 2 , where {\n)} is the 

computational basis. The expectation of the state over 



T&({2" 



-N/2 



\n)} n ) is given by 



E T (,,[^)(^r 2 ](x^|nn)(nfi| 

n 

^ \nifi) + \fhn) (nm| + (mn| 

= 2ng) m -^|nn)(nn|. 

This shows that the probabilistic mixture of 
T^* ) r ({2- Ar / 2 ,|n)}„) and a set of states {\n)} forms 
a complex projective 2-design. 

The resulting ensemble of states forms a complex pro- 
jective 2-design since the expectation of the state is pro- 
(2) 

portional to IIsym. Thus, the procedure to obtain a com- 
plex projective 2-design by using a phase-random circuit 
is given by 

1 . With probability 2 n +1 , choose a random bit n and 
generate a state |n). 



2. With probability 



2 N + 1 ' 



perform the CP phase- 



random circuit with an initial state !+ + •••+). 

This method of generating the complex projective 2- 
design requires O(N) random bits and 0(N 2 ) quantum 
gates in the CP phase-random circuit. The simple struc- 
ture of the CP phase-random circuit leads to an easier 
experimental implementation of the design than the pre- 
viously known results as listed below: 
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ancilla 



system 
bath 



l+> ■ 
l+> ■ 
l+> ■ 
l+> ■ 
l+> ■ 

l+> ■ 
l+> ■ 
l+> 
l+> 
l+> 



PRC 




QPE 



QPE 



FIG. 2: A quantum circuit to simulate canonical states. 
The upper half represents ancilla qubits and the lower half 
represents the system and the thermal bath. The PRC and 
QPE ( QPE -1 ) are abbreviations of the phase-random circuit 
and the (inverse of) quantum phase estimation, respectively. 



• In Ref . [J] , an exact unitary 2-design by using Clif- 
ford group is given, which requires 0(N a ) bits and 
0(N 2 ) quantum gates. In this method, we classi- 
cally calculate the unitary matrix and decompose 
it into one- and two-qubit unitary gates. Hence, at 
each time, it is necessary to repeat the calculation 
of the gate decomposition and to reconstruct the 
corresponding quantum circuit. 

• In Ref. Q, an e-approximate unitary 2-design is 
presented where the circuit is composed of one- 
and two-qubit Clifford gates and some gates are 
applied probabilistically. The number of gates 
is 0(N logl/e) in their definition of the 2-design 
corresponding to 0(N(N + logl/e)) in Defini- 
tion l[29|. 

• In Ref. [29j , an e-approximate unitary 2-dcsign by a 
random circuit [4J| composed of 0(N (N + log 1/e)) 
gates is proposed. All gates are randomly chosen 
from a so-called 2-copy gapped gate set, e.g., a set 
of controlled-NOT and single qubit rotations. 

• In Ref. [3l| . a local random circuit composed of 
0(Nt 4 (N + logl/e)) gates is shown to form an e- 
approximatc unitary ^-designs. The circuit is com- 
posed of random SU(4) gates acting on nearest 
neighbor qubits. 

In comparison with these results, the main advantage 
of our method for generating a complex projective 2- 
design is that it uses only diagonal gates and all gates 
can be applied simultaneously. 



2. 



Quantum simulator of canonical states 



The origin of the thermal equilibration of the system 
has been studied in the context of random states 
[23j | . Consider a system attached with a thermal bath 
and a Hamiltonian acting on the whole system, H to t = 
Hs + Hb + H int where Hs, Hb and H int denote the 
Hamiltonian on the system, that on the thermal bath 
and an interaction Hamiltonian between them, respec- 
tively. Let He be the Hilbcrt space restricted by a total 



energy e. It has been shown that, for almost any states 
of random states in H e , a reduced state in a system is 
close to the canonical state pth = e~ /3ffs /Tre~' 3 ' tfs where 
j3 is an inverse temperature determined by the total en- 
ergy e [20l | . It is also shown that this is the case for 
e-approximate complex projective i-designs [ioj . Such a 
typical equilibration is observed even for phase-random 
states although it depends on {r„, | e n )| whether the state 
equilibrates to the canonical state [22l l32l|. For instance, 
in the case of Y p hasc({-y=^, |e^) }) where d e — dim"H e and 
{|e,)} are the eigenstates of H tot in rl e , the reduced state 
in the system is almost canonical with high probability. 

Besides the fundamental importance on statistical me- 
chanics, these results offer a possibility to simulate canon- 
ical states by using random and phase-random states in 
H e , without knowing the eigenstates of the Hamiltonian. 
Particularly, the phase-random circuit can efficiently sim- 
ulate canonical states for classical Hamiltonians such as 
Ising models and classical spin glass models. Since all 
eigenstates of such Hamiltonians are separable, it is suf- 
ficient to prepare phase-random states in the computa- 
tional basis, which is approximately done by the phase- 
random circuit. 

To see this more clearly, let us consider N = N$ + Nb 
qubits where Ns and Nb are the number of qubits in the 
system and that in the thermal bath, respectively. Let 
Htot = $3* \hi){ni\ be a classical Hamiltonian where 
{|ni)} is the computational basis. We first define the 
restricted Hilbert space H e = span{|rij) |e — NS < Ei < 
e} where S is a small number independent of N. The 
total energy e determines the temperature of the system. 
For obtaining a canonical state, it is sufficient to generate 
TT p haso({-^=, which is approximately achieved by 

the quantum circuit presented in Fig. [5J 

The circuit is composed of r ancilla qubits initially 
prepared in the phase-random circuit, the quan- 

tum phase estimation (QPE) t 45i 4(3] and the projective 
measurement on the ancilla qubits in the computational 
basis. The qubits of the system and the bath are also 
prepared in |+)® , which is the equal-amplitude super- 
position of all eigenbasis of H to t- At each step, the state 
changes as follows: 

1. The phase-random circuit adds random phases 
{(fik} and generates a state in Y^({2 -JV / 2 , |rij)}). 

2. After QPE, the state is approximately |\& e ) = 



2 -n/2 j2 e 



l<Pk 



\Ek) where E k is a binary rep- 



resentation of the cigenencrgy E k . 

3. By performing the projective measurement 
P := {P e7 P-, e }, where P e is a projection 
onto H e and P-, e = I — P e , on the ancilla 
qubits, the state is probabilistically changed into 



\E k ) 



outcome corresponding to P e 
simulation fails. 



if we obtain the 
Otherwise, the 



4. The inverse of QPE changes the state into 



7 



(^T Efi t gw e e ^ fc ® 1+)*") so tnat we obtain a 
state contained in Y^({^i=, |nf)}) by tracing out 
the ancilla qubits. 

The error of the simulation mainly originates from 
QPE for two reasons. First, the eigenenergy is approx- 
imated by binary numbers within a precision of 2~ r . 
This approximation results in a round-off error j47| . 
Second, QPE does not exactly transform the state to 
l^e). This is inherited in the final state, resulting in 
— e) \^ f) + yfe \^ Error)- However, these errors can 
be suppressed by preparing a large number of ancilla 
qubits such that 2~ r <C AE where AE is the minimum 
energy gap of H to t- For local Hamiltonians, AE gener- 
ally scales exponentially with the system size TV, so that 
r = poly(N) is sufficient. 

Note that obtaining canonical states at low tempera- 
tures is generally difficult since the probability to obtain 
P e in the projective measurement depends on the tem- 
perature as shown in Ref. [23| . 

Our method is similar to the method proposed in 
Ref. [23| in the sense that both simulate the principle 
of equal weight on quantum circuits, from which canon- 
ical states are derived in statistical mechanics, and ex- 
ploit the projective measurement to obtain states in % e . 
Contrary to our circuit using only pure states, the cir- 
cuit in Ref. [23[ exploits mixed states and it works for 
any Hamiltonians. Comparing to the other methods for 
simulating canonical states [13, EH , our method may be 
simpler for experiments since all gates except the Fourier 
transformation in the phase estimation are diagonal in 
the computational basis and can be applied simultane- 
ously Thus, the experimental realization of the Fourier 
transformation immediately leads to the implementation 
of our circuit. However, our method has a disadvantage 
that it can efficiently simulate canonical states only of 
classical Hamiltonians. 



D. Sketch of the proofs 

In order to prove Theorem [1] and Theorem [21 we anal- 
ize how two-qubit diagonal gates in the phase-random 
circuits transform an initial state \4>o)- We denote the 
state after applying T two-qubit diagonal gates by \4>t) = 
Ut |0o)- We expand the state \<j>r) in the Pauli base 
and investigate the evolution of each coefficient. We de- 
note by p and q vectors corresponding to the subscripts 
of the Pauli base of two TV-qubit systems (pi, • • • ,Pn) 
and (<?!,••• , <?a0, where p it q t e {0, x, y, z}, respectively. 
Then, the state \<j>r) is expressed by 



\4>t){4>7 



r^?T(p,qK 



'q, 



where a p := a pi (g> • • ■ ® a PN and cr q .— Vqi *y - - - v qN 
are tensor products of the Pauli operators. Similarly, we 



(2) 

of a diagonal- unitary 2-dcsign , and expand \4> v ) in 
the Pauli base 

l</v)(0X 2 = z~ N J2 Up> q) CT P ® *q- 

p,q 

In terms of the expansion coefficients £r(p, q) and 
^(p,q), our goal is to show that for any initial state 



V(p,q) 



E Cr [fr(p,q)]- 



Mp,qX 



<22lv> ( 2 ) 



after sufficiently large T. By applying the similar argu- 
ment presented in Ref. [2^| , we can show that if Eq. ([2]) 
holds for any initial state, the phase-random circuit is an 
e-approximate diagonal-unitary 2-design. 

In Section IIV1 we present that the CP phase-random 
circuit achieves Eq. ([2]) for e = after applying 0(N 2 ) 
two-qubit diagonal gates, which gives the proof of The- 
orem [TJ In Section El we show that Eq. © holds for 
the CZ phase-random circuit after applying 0(N 2 (N + 
logl/e)) two-qubit diagonal gates, implying Theorem [2] 



IV. CP PHASE-RANDOM CIRCUIT 

We present a proof of Theorem [TJ To prove Theo- 
rem [JJ we follow the transformation of the expectation of 
an expansion coefficient E[£t+i (p, q)] by the CP phase- 
random circuit. 



A. Notations 

We introduce several notations that simplify our in- 
vestigations. Since the way how £t(p, q) is transformed 
depends on o~ p and <7 q , it is convenient to define sub- 
sets in {1, • • • , N} that specify the locations of a w (w = 
0, x, y, z) in a p and er q ; 



r (+) (p,q) 
r H ( P ,q) 

A(+)(p,q) 
A(->(p,q) 



= {* €{].,■•■ ,N}\ Pi = q t G 
= {«6{1,- - ,N}\ Pi = Qi G {x,y}}, 
= {%€{1,-- - ,N}]pi = qi€{0,z}}, 
= {»€{!,••• ,N}\ Pi = qie{0,z}}, 



consider a state 



Um \4>o) 1 where U v is an element 



where the bar sign of q~i represents to take a self-inverse 
'flip' map defined by = z and x = y. We denote the 
number of elements in each subset by the corresponding 
small letters, e.g., 7 < ' ± * > (P;q) (-^^(Pi q)) is the number 
of elements in T^^q) (A^^p, q)). We also denote 
the union of T (+) (p,q) and r (_) (p, q) (A (+) (p, q) and 
A (p, q)) by r(p,q) (A(p, q)). Similarly, the number 
of elements in T(p,q) (A(p,q)) is denoted by 7(p, q) 
(A(p,q)). By definition, 7(p,q) = 7 (+) (P, q)+7 (_) (P' q) 
and A(p, q) = A' + )(p, q) + '(p, q). For simplicity, wc 
often omit the part (p, q) in equations when there is no 
ambiguity. 
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We also define a function fs(p) of p, where S = 
{ii, ••• ,« s } is a subset of {1, • • • ,N}, such as fs(p) = 
(pi, • • • • • • ,Pi,, ' ' ■ ,Pn)- That is. the function f s 

flips all elements of p in the set S. For instance, 
fi,3(y,x, 0) = (y, a;,0) = (a;,a;,z). We also denote 
(/s(p),/s(q)) simply by /s(p,q). 

B. Calculation of E M (2) [^(p, q)] 

diag 

We calculate E w ( 2 > [£ v (p, q)] for an initial state \4>o) = 
£> n e iw » \fi) (Wn,r n > 0,ui n e [0,2tt) ). Since 

E ,y(2) K V (P,q)] = 2 _Ar E (2) [(^IfTp |^) (^IfTq |^)], 

"diag M diag 

some calculations lead to the coefficients. For (p, q) sat- 
isfying 7(p, q) + A(p, q) = TV, we obtain 

E W ( ? >jap,q)] = 2- N £« H s njmj 

x 5 nm (1 - 2rij) 

x n (i-^™,) 

x ]J -*( 1 -^m J )( <5 P^- (S P J i;)( 1 - 2n j)' ( 3 ) 

where S nm is Kronccker delta and nxn^-'-n^ 
(rrixTrii ■ ■ ■ ton) is a binary representation of n — 1 



(to — 1). In other cases, E < 2 ) [£ v (p, q)] = 0. 



C. Transformation of E c cp [£r(p, q)] 

We consider the transformation of a state by the CP 
phase-random circuit. Let us consider the expansion of 
E c cp[|0t)(0t|® 2 ] in the Pauli base given by 

E C?P UtH^tH = E E c ? r [6r(p, q)]<7 p ® er q . 

pq 

For simplicity, hereafter we omit the subscript 
for the expectation values. By applying 
Wt+i {q.t+1 , Pt+i ! 7t+i) on a pair of qubits specified by 
two numbers (i, j), the expectation of coefficients changes 
to 

E[£ T +i(p,q)] =2~ 2W £ G y (p,q;p',q')IEKT + i(p',q')], 

(P',q') 

where the matrix Gy(p, q; p', q') is given by 

Gij{p, q; p', q') = E[Fra 9 WT+iap>w£ +1 Tr<rJVT+i(r qL >W}. +1 ]. 

In Appendix [XJ we present a derivation of the matrix 
Gij(p, q; p', q')- Then, we obtain 



E[fr+i(p,q)] = { 



E[£r(p,q)] 

i(E[eT(p,q)]±E[e T (/i(p,q))] 

+EKr(/ i (p > q))]±EK r (/y(p > cO)]) 

±(E[fr(p,q)] + vE[Z T (f j (p, q))] 

+ U E[£ T (/ I (p,q))]+m;E[£ T (/ y (p,q))]) 



if i,j £ A(p,q), 

if J er( ± )(p,q), J e A( P ,q), 
ifier(")(p,q),ierM(p,q), 

otherwise, 



(4) 



where u and u is ±1. 

Note that a set of A(p, q) and a set of T(p, q) are both 
invariant under the transformations since the transfor- 
mations are composed of a function /s(p, q) which flips 
x (y) to y (x) and (z) to z (0). 

From the last case of Eq. (QJ, for any pairs of in- 
dices (p,q) satisfying 7(p,q) + A(p,q) < N, we have 
E[£r(p, q)] = for any T once after i £ A(p, q) U T(p, q) 
is chosen. The other cases show that when one of 
is in r(p,q), we take the uniform average of the 
'flipped' terms. In particular, when i € T^^p, q) is 
selected, the term acquire a negative sign. Hence, af- 
ter all pairs of qubits are drawn, the expectation of the 
state becomes an average of all flipped terms with ap- 



propriate negative signs. To express the expectation, 
it is convenient to introduce ri^J„(p,q) and r^(p, q) 
as subsets of r(~)(p,q) of which the number of ele- 
ments is even and odd, respectively We also intro- 
duce S even (p, q) and ^^(p, q), which are defined by 
5 el , e „(p,q) = {/ s (p,q) where s C A(p,q)Ul ,+ )(p,q) U 
rLe„(p,q)} and 5 dd(p,q) = {f s (p, q) where s C 
A(p, q) U r< + ' (p, q) U r^](p, q)}. Then, we obtain the 
following proposition. 

Proposition 2 After applying Wij to all combinations 
of i and j, which requires Tap = jV ^~ 1 ^ two-qubit gates, 
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the coefficient converges to 
E[6r r(p,q)]=2-" 



E E 

S c „cn(p,q) Sodd(p,q) 



£o(p',q'), 



(5) 

where the summation is taken over all (p',q') G 

S'e'uen(odd)(P7q)- 

Finally, by calculating Eq. explicitly, we obtain that 
for any initial states and for all (p,q), 

E[£T CP (p,q)] =E (2i [C v (p,q)]. 

M diag 

For the detailed calculation, see Appendix |B] This con- 
cludes the proof of Theorem [TJ 



V. CZ PHASE-RANDOM CIRCUIT 



The second lemma states that the convergence time 
Tconvit) defined in Theorem [5] scales cubic of the system 
size N . 



Lemma 2 For any initial state, the convergence time 
T C onv(f) satisfies 

A 3 N 2 

— \og2 + —\oge- 1 +0(N 2 )<T conv {e) 

< 3A 3 log 2 + A 2 log <T X + 0(N 2 ). 

Therefore, T conv {e) = 0{N 3 ) + 0(N 2 ) loge" 1 . 



We present the proof of Lemma [T] in Subsection IV Al and 
the proof of Lemma [2] in Subsection IV Bl 



Theorem [5] is shown by proving the following two lem- 
mas. The first lemma guarantees that the state trans- 
formed by the CZ phase-random circuit for any initial 
state converges to the corresponding toric 2-design. 

Lemma 1 For any initial state and V(p,q), 



lim 

T->oo 



E c? z [£t(p, q)] - E zy( 2 > Mp> q)] 



A. Convergence of the distribution 

1. Transformation of E c cz [£t(p, q)] 

We consider the transformation of a state by the CZ 
phase-random circuit. Similarly to the CP case, we ob- 
tain 



E[£ T +i(p,q)] 



(E[£ T (p,q)] A(p,q), 

KEIM/ifaq))] ±E[fT(/< 3 -(p,q))]) if i G r(±)(p,q), j G A(p,q), 
J(E[fr(p ) q)]+«E[f r (/i(p,q))] 

+uE[€r(/i(p,q))] +««E[fT(/«(p,q))]) if i G r(")(p,q), j G rW(p,q), 

otherwise. 



(6) 



where u,v E {+, — } (sec Appendix |A"|) . The transfor- 
mation is different from that of the CP phase-random 
circuits given by Eq. Q only when i G r^^q) and 
j G A(p, q). This prevents the circuit from efficiently 
randomizing the corresponding phases. Thus, we have 
to introduce stochastic transformations by choosing (i,j) 
randomly for achieving a diagonal-unitary 2-dcsign. 



2. Modified phase-random circuit 

Similarly to the CP case, for any pairs of indices (p, q) 
satisfying 7(p, q) + A(p, q) < A, we have E[£t(p, q)] = 
for any T once after i ^ A(p, q)UL(p, q) is chosen. In or- 
der to avoid the complication by dealing with such (p, q), 
we first apply a two-qubit gate W^j (a, ft) on all neighbor- 
ing qubits (2k — 1, 2k). We denote this unitary operations 



by W = Wn,n-iWn-2,n~3 ■ ■ • W 2 ,i- (see Fig. [3]) When 
A is odd, we define W by Wn,n-iWn-x,n-2 • • • W^,!- 
The number of the two-qubit gates required to perform 
W is Ty^ = [A/2] where \n\ is the smallest integer larger 
or equal to n. Note that W is composed of commutable 
gates and is deterministic. Hence, they can be applied 
simultaneously. 

In an analogous way with Proposition^ we obtain the 
following proposition. 

Proposition 3 If 7 (p, q) + A(p, q) < A, E[£ T ^ (p, q)] = 
0. If (p, q) satisfies 7(p, q) + A(p, q) = A, 



E[&v(p,q)] 



E E E m 



■.crfci, 



scr 



s > ° /A(p>q)), 



«'cr(+) 
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r(h,jl) TT/(i2J2) 









w 















FIG. 3: A modified phase-random circuit. A iV-qubit 
unitary operation W is applied in advance of the phase- 
random circuit. The unitary operation W is composed of 
two-qubits diagonal gates acting on neighboring qubits, W = 

Wjv,JV-1 ■ ' ■ W3,2W2,1. 



where (p, q) is omitted for simplicity and A(p, q) is a set 
of i G A(p, q) such that i is paired with an element of 
r(p,q) in W. 



3. Transformation o/E[£t(p,q)] after applying W 

By using the modified phase-random circuit, after ap- 
plying W that consists of two-qubit gates, addi- 
tional T' two-qubit diagonal gates randomly selected 
from Wdiiig are applied on a randomly chosen pair of 
qubits. The following proposition shows how E[£r(p, q)] 
relates to E[£ T+1 (p, q)] for T = T W + T' where V > 0. 



Proposition 4 For T > T, 



E[6r +l (p,q)]= J2 £(p,q;p',q')E[Mp',q% (7) 

(P',q') 

where £(p, q; p', q') is equal to A(A ^^i]"^ if (p'i q') = 

(p, q), n (n-i) if (p'» i') = /*(p> i) for * e A (P' <*) and 

otherwise. 



Proof 2 The equation ([7]) is equal to 

E[fr + l(p,q)] = A(A " 1 ^ 7 ir 1) ^T(p,q)] + 



N(N-1) 



2l 



N(N - 1) 



^ EK T (/ 4 (p,q))]. (8) 

eA(p,q) 



We show Eq. © recursively. 

Firstly, wc show that E[£t- +i(p, q)] = 
E [&vl/i(p,q))] when i G A(p, q) and j € r(p,q), 
and E[£r^+i(p,q)] = E[£r^(p,q)] otherwise. This 
is derived by using Eq. and Proposition [3] in the 
following way. When i,j € A(p, q), Eq. © implies 
E[fr-+i(p,q)] = E[e T , v (p,q)]. When i e A(p, q) and 



j G r^^p, q), E[^7^+i(p, q)] is calculated to be 
EK7V+i(p,q)] 



^E[CT ri ,(/,(p,q))]±iE[eT,,(/ y (p,q))] 



= E[^(/ i (p,q))] J 

where the second line of the equation is obtained from a 
property of E[£x^, (fi{p, q))] such that 

Vj G rW(p, q), E[Zt» (fj(p, q))] = ±E[^ (P, q)]. (9) 

When i,j G T(p, q), a direct calculation shows that 
E[&-+i(p,q)] = E[^(p,q)] by using Eq. ©. Thus, 
we obtain the statement. 

Since the probability that j G T(p, q) is chosen for 
a fixed i G A(p, q) is given by 7/(2)1 where („) is a 
binomial coefficient, Eq. © is shown for T = T-^. More- 
over, E[^.+i(p,q)] also satisfies E[^+i(/*(p,q))] = 
±E[f T +i(p,q)] for » G r( ± )(p,q), so that Eq. © is 
recursively obtained. 



Proposition S] implies that, for T > T^, E[£r(P)Q)] is 
given by a convex sum of E[£r^.(/ a (p,q))] where s is a 
subset of A (p, q). Hence, if wc specify L^' C {1, ■ ■ • ,N} 
where n L^' = 0, the transformation of E[£y(p, q)] 
is closed in S(l' + ',l("') := {(p, q)|A(±) (p, q) = 
L^)}. In Proposition [5J we consider the transforma- 
tion in Y,(L( + > , and derive the stationary distri- 
bution E[^ QO (p,q)] := lim T _ +00 E[^ T (p, q)] for (p, q) G 
£(£(+), £(")). 

In order to obtain Proposition [51 we use the Perron- 
Frobenius theorem (49j for irreducible and aperiodic non- 
negative matrices M . Irreducibility is a property such 
that, for all i and j there exists a natural number n 
such that (M n )ij > and aperiodicity is a property that 
Mjj > for alH. A non- negative matrix implies that for 
Mij > for all i and j. The Perron- Frobenius theorem 
is given by the following statement. 

Theorem 3 (Perron- Frobenius theorem [49]) If a 

non-negative matrix M is irreducible and aperiodic, the 
maximum eigenvalue A > is uniquely determined. Let 
j A) be the corresponding eigenvector to the maximum 
eigenvalue, then, lim„_ J . 00 (iM) n = |A) (A|. 

In addition to the irreducibility and the aperiodicity, 
when a non-negative matrix M is bistochastic, that is, 
'^2, i Mij = y"V Mij = 1, it is known that the maximum 
eigenvalue A is equal to 1. By applying these facts, wc 
obtain Proposition [3] 

Proposition 5 Let L^> be a proper subset of 
{1, •••,#} satisfying n = and L := 

IWULH £ {1,---,N}. For (p,q) G EfL'+U 1 "'), 
the stationary distribution E[^ 00 (p, q)] is uniform in 
E(£(+), £(-)), that is, 



■Koo(p,q)] = 



53 E[^(p',q')], (10) 

(p',q')eS(L(+),i(-)) 
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where I is the number of elements of L which is also equal 
to A(p, q). Moreover, for any (p, q), we obtain 



E[£co(p,q)] =E ( 2, Mp,q)] 

diag 



(11) 



Proof 3 In the case of L = {1, • • • , N}, it is straight- 
forward to show Eq. (|TTj) from Eq. ([5]). When 
L 7^ {l,--- ,N}, we can obtain Eq. (fT0|) by ap- 
plying the Perron-Frobcnius theorem to the matrix 
Q in Y,(L( + \ L^~~>). If we restrict the matrix Q in 
£(l/ + ), L/('), it is straightforward to see that Q is an ir- 
reducible, aperiodic and bistochastic non-negative matrix 
in L('). Hence, the Perron-Frobenius theorem 

guarantees that there exists a unique stationary distribu- 
tion in £(X(+), £K>). Since the evolution governed by Q 
in £(i/ + ) , L^') is uniform, the stationary distribution 
is also uniform, leading that, V(p,q) € L^), 



-Koo(p,q)] = ^7 



£ E^p'.q')]. 

(p',q')e£(L(+),L(->) 



Thus, we obtain Eq. ([TO]). 

Next, we show Eq. (fTTj) . For (p, q) satisfying A(p, q) + 
7(p, q) < N, E^oc (p, q)] = since all E[£ T ^ (p, q)] in the 
right hand side of Eq. (fT0|) are zero as shown in Proposi- 
tion [31 Since Eq. implies E[£ u m (p, q)] = for such 

diag 

(p, q), we obtain E[£ Too (p, q)] = E[£ (2) (p, q)]. 

diag 

When (p, q) satisfies A(p,q) + 7(p,q) = N, we sub- 
stitute E[^t^,- (p, q)] given in Proposition [3] into Eq. (fTUfl . 



and obtain 



!K<x,(p,q)] =2 



— N 



E E 

S c „ e „(p,q) S 0(Ili (p,q) 



£o(p',q')- 



As shown in Appendix IBl this is equal to E^pj [£ v (p, q)] 



B. Convergence time for the phase-random circuits 

In this section, we investigate the convergence time 
Tconv(e) defined by VT > T conv (e), 



l c cz[U® 2 ®U\ 



, P j [L/| 2 ®[/t» 2 ]| 0<e . (12) 



As explained in Section IIIID1 a sufheient condition for 
Eq. (fl2|) to hold is given by 



V(p,q), 



E c cz Kt(p, q)] -E c c Z Koo(p, q)" 



< 2537- (13) 



Note that E c cz (p, q)] = E (2 ) [£ v (p,q)] from 

°° diag 

Lemma [TJ Similarly, we can obtain a necessary condi- 
tion for Eq. (fT2"j) to be satisfied by evaluating a lower 
bound of the diamond norm: 



V(p,q) 



E c£ z (p, q)] - E c cz [Coo (p, q)] 



< e. 



(14) 



We derive an upper and a lower bound of T conv (e) 
by using Eqs. (fl~3|) and ([T4| . respectively, and prove the 
Lemma [2] stating that for any initial state, 

/V 3 N 2 

— \og2 + —\oge- 1 +0{N 2 )<T conv {e) 

< 3N 3 log 2 + iV 2 log e- 1 + 0(N 2 ). 

Note that, for (p, q) satisfying A(p, q) + 7(p, q) < N, 
T co „v(e) <T W = \N/2\ since E[£r(p,q)] = for T > 
Tyy, and, for (p, q) satisfying A(p,q) = N, T conv (e) = 
since E[£r(p, q)] = £,o(p, q)- In the following, we consider 
only (p, q) such that 7(p, q) + A(p, q) = N and A(p, q) ^ 
N. 

In order to show Lemma [21 we use techniques of the 
Markov chains (50[. We provide a brief introduction of 
the Markov chains in Appendix [C] We map the trans- 
formation of E[£r(p, q)] into a Markov chain and give a 
lower and an upper bounds of the convergence time. 



1. Map to a Markov chain 

We present a map from the transformation of 
E[£oo(p, q)] by the (modified) CZ phase-random circuit 
to a Markov chain. Since the transformation in the 
phase-random circuit is deterministic for T < T^, we 
consider the evolution for T > as a Markov chain. As 
shown in Proposition^ the transformation of E[£x(Pi q)] 
is closed in L^) and Q [(p, q); (p', q')] satisfies 

a Markovian property. Moreover, it is observed from 
Proposition S] that G[(p, q); (p', q')] i s equivalent to a 
transition matrix of a random walk on a /-dimensional 
hypcrcubc where each vertex is given by (p, q) G 
E(L( + ), L^ - '). Note that Q is irreducible and aperiodic in 
T,(L^ + ',L^~'). However, E[£r(p, q)] cannot be regarded 
as a probability distribution of a Markov chain since they 
are not necessarily to be non-negative. Instead, we define 
a probability distribution in the following way. 

We set the initial probability distribution {n (p,q)} 
on the Markov chain M{L^ + \ L^) to be 



n (p,q) 



E[gr tV (p,q)]-n min (L(+>,L(->) 
n sum (£(+),£(")) 



where 



n min (L«, min E[^( P ,q)], 



and 



n sum (Z« L(-)) := 

J2 (E[^(p,q)]-n min (L(+),L(-)) 



When there is no ambiguity, we omit the variable descrip- 
tion (L( + \ L(') for II m j n and II sum . Then the probabil- 
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(0,0,0) 



(0, 


3,1) 




(1,0,0) 


y 




(0,1,1) 





(1,0,1) 



(1,1,1) 



(0,1,0) 



(1,1,0) 



FIG. 4: A 3-dimensional cube. A random walk on the cube 
is equivalent to a Markov chain M with N = 3. 



ity distribution {IIi(p,q)} is calculated to be 
ni( P ,q)= Yl ^[(p,q);(p',q')]no(p',q') 

(p',q')6£(L( + ),L<->) 
1 



n s 



-(E[^+i(p,q)]-n roin ), 



where we use Proposition|4]and a fact that the matrix Q is 
bistochastic. Repeating this, the probability distribution 
after k steps is obtained by 



n fe ( P ,q) 



1 



-(E[^+ fe (p,q)]-n min ). 



(15) 



Thus we can define a Markov chain A4(L^ + \ L^) 
on a Z-dimensional hypercube with a transition matrix 
G[(p, q); (p', q')] an£ l a probability distribution life (p, q). 
Note that Eq. (HU) leads to V(p,q) G £(")) 



nfc(p,q)-noo(p,q) 



E[e T ^+fc(p J q)]-EK 0o (p,q)] 



where n oo (p, q) is a stationary distribution of the 
Markov chain 7W(L (+) , L (_) ). This implies that if the 
Markov chain M(X (+) , converges with an error 

e/n sum , E[^t^(p, q)] converges with an error e. Hence 
Tconv(e) = T mix (e/n sum ) where T mix is the mixing time 
of the Markov chain defined in Appendix. [C] 

The mixing time of the Markov chain on a hypercube 
depends on two factors. One is the dimension of the hy- 
percube I and another is the probability that no change 
happens on a Markov chain, which is referred to as a 
staying probability. Obviously, a larger I and a smaller 
staying probability result in the longer mixing time. In 
the Markov chain Ai(L^ + \ L^'), the maximum of the 
dimension and the minimum of the staying probabil- 
ity are achieved for I = N — 1. If the Markov chain 
A4(L( + \ L(~>) with I = N — 1 converges with an error 
e/n sum , the other Markov chains with I ^ N — 1 con- 
verges with an error less than e/II sum . Thus, hereafter, 
we consider only the Markov chain M.{L^ + \ L^) with 
I = N — 1 , which we denote by Ai . 



Finally, we simplify the notations of the Markov chain 
M.. Since it is equivalent to a Markov chain on a (N— 1)- 
dimensional hypercube with a transition matrix Q, we la- 
bel each vertex by a binary number i = ■ ■ ■ , in-i) G 
{0, l}-^ -1 , not by (p,q). If we identify (po,qo) with 
io = • • • 0, a new label for /fe(po,qo) is given by 
■ • • 010 • ■ • where only the fc-th digit of the binary rep- 
resentation is 1 . The vertices i and j are connected if and 
only if H(i, j) = 1 where H(i,j) =J2\ik~jk\ (sec Fig. [p. 
In this new labeling, the transition matrix Q[p, q); p', q')] 
is simplified to 



i 



N(N- 





if i = j, 
jy HH(i,j) = l, 

otherwise. 



2. 



Lower bound of the mixing time of a Markov chain M 



The Markov chain on a {N — l)-dimensional hyper- 
cube with a staying probability 1/2 has been well stud- 
ied. For such a Markov chain, the transition matrix 
is given by 1/2 + V' /2 where V is a matrix of which 
elements are jy±t f° r = 1, and otherwise. 

All eigenvalues of the transition matrix are known to 
De {1 — iv^i }fc=o,--- ,n-i H^. On the other hand, it 
is observed from Eq. (HU) that V(i,j) = (1 - + 
jrV. Hence, the eigenvalues of V{i,j) are given by 
i 1- jv(jv-i) } fc =o,---,jv-i, which leads to max i=2 ,3,- |A»| 
1 



N(N-l) JK=u,-,JV-l, 

wj^irn • Using Eq. (|C1|) in Appendix [Cl the mixing 
time T m i x (e) is bounded from below by 



,N(N — 1) . 1 

(-^ --l)log^<T mix (e). 



3. Upper bound of the mixing time of a Markov chain M 

In order to derive an upper bound of the mixing time, 
we slightly modify the transition matrix by chang- 

ing the the staying probability from 1 — to 1 

A new transition matrix V(i,j) is given by 



N-l ■ 



1 



1 

N-l 







if i = 3, 
i£H(iJ) = l, 

otherwise. 



We denote by M. a new Markov chain with a transition 
matrix V . Since the staying probability of V is smaller 
than V , the mixing time of M provides an upper bound 
of that of M . 

We investigate the mixing time of the Markov chain 
SA by the coupling method (see Appendix [C]) . For con- 
structing a coupling, we interpret M. as follows. At a step 
t = k + m(N-l) {k = !,■■■ ,N-1 andm = 0, ■), 
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choose j G {ii, ■ • • , in-i} at random and, if j = i ko , flip 
ik , otherwise, do nothing. 

Based on this, we define a coupling (X k , Y k ). A state in 
Xk (Yk) is a bit sequence x\ ■ ■ ■ xn^i (yi ■ ■ ■ Vn-i) where 
Xj (y/) e {0,1}. At a step t = k Q + m(N - 1) (fc = 
1, • • • , N — 1 and m = 0, 1, • • • ), we randomly choose 
j G {ii, • • • ,ijv-i}- If j Hoi do nothing. If j = i fco and 
= yj , flip both of them. If j = i ko and a;-,- 7^ ?/j , flip 
Xj and where a(j) is the position satisfying x a u\ 7^ 
2/ a (j) next to j. For instance, when X = (0, 1, 0, 1, 1) and 
Y = (1,1,0,0,0), a(l) = 4 and a(4) = 5. Each of X fc 
and Yk is equivalent to the Markov chain A4 and they 
also satisfy X k +i = Yfc+i if X k = Y k . Hence, (X k ,Y k ) is 
a coupling of M. 

For investigating the coupling, we use the property of 
a special type of Markov chains, a coupon collecting. A 
coupon collecting of r-coupons is a Markov chain on a set 
of states {0, 1, • • ■ , r} with the transition matrix given by 

ProbLY t+1 = k + l\X t = k] = — -, 

r 

Pxob[X t+1 =k\X t = k} = -. 

r 

A coupon collecting of r coupons is interpreted as a trial 
of collecting a complete set of r different coupons by 
drawing one coupon at each step. When ro coupons are 
initially at hand, we denote the necessary number of steps 
to draw all coupons by T^oupum- 

Proposition 6 Let (X k , Y k ) be a coupling of M. defined 
above and T xy be its stopping time. Then T xy is bounded 
from above by 

T < (N- i) T (x-hN-2) 

J-xy Zi -Ly 'coupon 

Proof 4 By definition of the coupling, for a fixed ko £ 
{1, ■ ■ • , N — 1}, once i ka is chosen at the steps k = 
k + rn(N - 1) (m = 0,1,---), we have x ikQ = y lkg . 
The number of steps necessary for picking up ko from 
{1, • ■ ■ , AT — 1} is equal to that of a coupon collecting of 
(N — l)-coupons with initially N — 2 coupons at hand, 

that is, rio^pon^ 2 ^. Since there exists N — 1 choices 
of fco, the stopping time T xy is smaller than or equal to 

(N- lW (W " llW " 2) 

X 1 » ± ) ' coupon 

■ 

A coupon collecting of r-coupons starting with no 
coupon is a well-studied problem. It is known that if 
we draw coupons r log r times, we can complete all of 
r-coupons with a high probability, that is, rioupon is typ- 
ically given by rlogr. For rijuponj we first show that 

Prob^l > r(Iog(r - r„) + c)] < e" c , (16) 

for any c > 0. This is shown in a standard way (see 
[D]for details). By appealing to Theorem @] presented in 



Appendix [Cl Proposition [6] and Eq. (|16[) , we can bound 
A((iV-l) 2 c) by 

A((N - ifc) < maxProb^y > (N — lfc)\ 

x,y 

< Prob[(7V - l)ri™- 2 > > (AT - l) 2 c] 
= Prob[r(™- 2 )>(7V-l)c] 

<e" c . 

Denote by T m j x (e) the mixing time of the Markov 
chain M with an error e. Since T m j x (e) is defined by 
A(T m i X (e)) < e, we obtain 

f m Ue) < (AT-lfloge" 1 , 

which also provides an upper bound of the mixing time 
of the Markov chain Ai since T m j x (e) < T m i x (e). 

4- Upper and lower bounds of the convergence time T conv (e) 

We require an error e/(II sum 2 2Ar ) in order to obtain 
an upper bound of the convergence time T C onv(e) (see 
Eq. (I13p ) and, an error e/II sum for a lower bound (see 
Eq. (Tj"4")) h Recalling that the unitary operation W con- 
sists of Ty^ = \N/2] two-qubit gates, we obtain bounds 
for T conv (e) such that 

( 1 l)!og( Y e ) - T «>"v(e)-% 

<(N- l) 2 log( - suml '— j ). 

Although n sum (L( + ), L(~}) depends on an initial state 
\(j> ), it is shown that n sum (L ( +\ i ( ~)) < 2 N for any 
(L^ + \ L,(~>) and any |</>o). Therefore, we finally obtain 

/V 3 A^ 2 

— log 2 + — log e- 1 + 0(N 2 ) < T conv (e) 

< 3A 3 log 2 + A 2 log e- 1 + 0(A 2 ), 

for any initial state |0o)- This concludes the proof of 
Theorem [2j 



VI. SUMMARY AND CONCLUDING 
REMARKS 

In this paper, we have introduced concepts of diagonal- 
unitary i-designs and toric f-designs that simulate up to 
the t-th order of statistical moments of diagonal-unitary 
matrices and phase-random states, respectively. We have 
presented how to implement diagonal-unitary 2-designs 
in the computational basis for iV-qubit systems by using 
two types of the phase-random circuits, the CP and the 
CZ phase-random circuit. We have shown that the CP 
phase-random circuit exactly achieves a diagonal-unitary 
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2-dcsign after applying 0(N 2 ) up to two-qubit diagonal 
gates. On the other hand, the CZ phase-random circuit 
approximately achieves a diagonal-unitary 2-dcsign after 
applying 0(N 2 (N + logl/e)) two-qubit diagonal gates, 
showing that random variables in the genuine two-qubit 
diagonal gate provide stronger ability for randomizing 
phases. Since they are composed of only diagonal gates 
and, for the CP phase-random circuit, all the gates can 
be applied simultaneously, it is expected that the phase- 
random circuits are practically easy to implement. 

We have also provided applications of the phase- 
random circuits. One is a generation of complex- 
projective 2-design by combining a phase-random cir- 
cuit with a simple classical procedure. Compared to 
the previously known results, our method is simple for 
implementation in the sense that all gates are diagonal. 
Another is a quantum simulator of canonical states of 
classical Hamiltonians where canonical states with any 
temperature can be implemented by preparing a proper 
initial state, without realizing the Hamiltonian itself. 

In analogy with the random circuits, which are shown 
to approximately achieve unitary ^-designs for any t by 
applying poly(N,t) two-qubit gates [3l|, it is natural to 



guess that the phase-random circuits with appropriate 
gate sets would also achieve diagonal-unitary i-designs 
in poly(N,t) iterations. However, it is not the case if we 
use the gate sets composed of only two-qubit diagonal 
gates since there is a lack of the number of parameters 
due to the commutability of gates. Thus, for construct- 
ing diagonal-unitary ^-designs, the gate set should in- 
clude multi-qubit gates if only diagonal gates are used. 
It would be interesting to specify the diagonal gate set 
of the phase-random circuit achieving diagonal-unitary 
i-designs, or to construct a quantum circuit composed 
of non-diagonal two-qubit gates which achieves diagonal- 
unitary i-designs. 
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Appendix A: Calculation of Gi 3 (p, q; p', q' 



Gjj (p, q; p', q') is defined by 



Gn (p, q; p', q') = E[R«r p Wy o^Wj I &o- q W«o- q 'Wj] 1 

where W%j is a two-qubit diagonal gate on the z-tli 
and j-th qubits and randomly chosen from the gate 
set W cz = {diag(\,e la ,e lf3 ,~e l ( a+ ^)} atl3 or W CP = 
{djoj(l,e la ,e^,e 17 )} ai) g j7 . Since W CP is more gen- 
eral than W cz , we start with the calculation of 
G«(p,q;p',q / ) for W CP . 

In order to calculate Gjj(p, q; p', q'), we define T> ao , £ ao 
and A Q by 



T^ab = $a0$bz + S az SbO, 
&ab &ax&by SaySbxi 

and 

A a = S a + S az — S ax — Say 



Acknowledgement 
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comments on unitary i-designs and P. S. Turner for hclp- 



Wc also use a notation that S n& s = 1 if n S S and 
<5„ eS = if n S. 

For Wij = diag(l, e ta , e 4 ^, e* 7 ), it is straightforward to 
calculate Trap Wiy-ay and we obtain 
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■^TrapWijdp'Wl = <w{<Wie{o,z} + 5<W{*,i/}(cos£ + cos(a - 7)) + ^S Pje{x , y} (cosa + cos(^ - 7)) 



2^Pi-Pj£{x-y} [ cosa + cos/3 — COS7 — cos(a — {$) + cos(a — 7) + cos(/3 — 7)] 



- [<W{o,z}(shia - sin(^ - 7)) + 6 Pi€ { Xty y (sin 7 + sin(a - 0))]£ Pj , v f. 

+ \ S Pj-p'j |^e{x,j/}( C0SQ! - cos (^ - l)) v Pi;P[ 

- fee{o, 2 }(sin^- sin(a - 7)) + <5p je {a,,y} (sin 7 + sin(/3 - a))]S Pi y. 
- -(sin/3 + sin(a - 7))^^;^^ - ^( sin « + sin(/3 - 7))^p<,p^Pj, P J 



(Al) 



r 



In the case of W c ' p , by taking the average over a, /3, 7 = 0, 4^, is calculated to 
I 



p G y (P'*P il) = <W'^q,q' (^Pi.Pi.gi.Qye^^ + ^pi.gie^ 

+ ^Pi,p'i^qi,q'i ^P3,9i6{o,z}^pj,p^« 3 ,g< ^p j + ^Pj ,p'j^qj ,q'j ^ < V;,9;e{o,z}^Pi,p-'5<?i,^ + ^PiiPj^g.,^^ 

4- i/ 7 /£ /£ z/ 7 

T 4°Pi.Pi°9i,^ Pj,Pj 9j,9j- 



By investigating each case, we obtain Eq. (HJ). 

On the other hand, in the case of W cz , 7 is set to be 



a + /3 + 7r and Cry is obtained as 



G ii (p,q;p / ! q / )=2 



2N 



^P.P'^q.q' {&Pi,Pj,qi,qje{0,z} + . ^P;,Pj,<Ji,?jG{a;,!/}) 



+ 2^ P " qi£ { X ' V }^ P " P 'i^ q " ql i^ P i' P 'i^ >q j' q 'i 2^ Pj ' P 'i^ qj ' q '^ ^ ?^Pj' q j G { X 'y}^ P j> P 'j^ q J' q 'j(^ P i> P 'i^~ )q i> q i ^ 2^ P *' P *^ 9i,9 i' 

-I T) ,73 ,£ z/ 7 , -4 Z 7 'Z 7 ,73 ,73 / 4- i/ 7 ,£ if if 1 

' 2 Pi'Pi H' q i Pi >Pj q 3 'Ij ^ 2 Pi,P i qi ' q i Pi'Pj Ullj 1 ^PUPi<hAiPj,Pj<lj,<lj 



r 



leading to Eq. 



Appendix B: Calculation of the expectation of the 
state after the phase-random circuits 



Here, we show that 



-N 



E E 

'S e -»en(P,q) S odd (p,q) 



Co(p',q') = E w (2, [^(p,q) 



(Bl) 
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For an initial state of the phase-random circuit rep- 
resented by |0 O > = E r « el ""l">> £o(P,q) is given by 
£o(p,q) = 2 _JV (0 O | cr p |0o> (0okq|0o>- A straightfor- 
ward calculation shows that 



^(P^q)^ 2 N E r ^ r m r irke % 



i{w n -w m +ui—w k ) 

n' m' l< k ' 

x S njmj Si ]kj [5 Pj0 + S PjZ (l - 2rij)(l - 2lj)] 

j€A(+) 

x S njmj Si ]k] [S Pj0 (l - 2lj) + 8 PjZ (l - 2rij)} 

x Y[ 3n ]m] 5i ]k] [S PjX - S PjV (l - 2rij)(l - 2lj)] 

x n i{8 nimj 5 hkj [S Pj x(l - 2l i) + W 1 - 2n i)l- 

jer(-) 

where 8 nm = 1 — <5 nm . By substituting £o(p?q) into the 
left-hand side of Eq. (|B1[) , we obtain 

EKoo(p,q)] =2- 2JV ^ W,^-""^-"') 
x 28 njmj 8i j k j 8 nj i j 28 njmj 8i j k j S nj i j {l - 2rij) 

jer(+) 

x -2i8 njkj 8 mj i j {8 PjX - <5p 3 y)^„ imj (l - 2n 3 -), 

jer(-) 

where we have used a relationship that 



E - E 

Set>en(p,q) S odd (p,q)' 

jer(-) 

= JJ -2iS njmj 8i jkj (8 PjX - 8 Pj y)(nj - 

= -2i8 n]kj 8 m]l] {8 P]X - 8 Pjy )8 njmj (1 - 2nj). 



Finally, using A<+> + A (_) +7 (+) +7 (_) = AT, we obtain 
the result such that 

2-4 E - E W,q') = 

S , e U „n(P,q) S ocid (p,q) 

2- w E r «^ n ^ n ^^-(1-2%) 

n,-m jeA<+) j'eA(-) 
x 8 nm -i8 nj m j (Sp :j x — 8p j y)(l-2nj), 

which is equal to E ( 2 > [£ v (p, q)] given by Eq. 

diag 



Appendix C: Introduction of Markov chain 

A Markov chain is a sequence of random variables in- 
dexed by a discrete step t £ N that take values in a 
set of states S = {s}. We define a probability distri- 
bution {IIi(s)} se s at a step t over the state space 5*. 
The Markov property is that a probability distribution 
n t+ i depends only on IT. This evolution of the proba- 
bility distribution is governed by a stochastic transition 
matrix V such that Ht+i = "PIT. The element of a tran- 
sition matrix is denoted by V(s,s'), which represents a 
probability that a transition from ,s to ,s' occurs. Using 
an initial distribution IIo, the probability distribution at 
step t is given by ±l t = 'P'llo. A Markov chain is said 
to be irreducible (aperiodic) when a transition matrix 
is irreducible (aperiodic). For an irreducible and aperi- 
odic Markov chain, the Perron-Frobenius theorem guar- 
antees that there exists a unique stationary distribution 
ILxj = limt-j.oo II t independent of the initial probability 
distribution. 

We define the mixing time. The mixing time is the 
number of steps required for the actual distribution to be 
close to the stationary distribution, where the distance 
after i-steps is defined by 

A(t) :=max|n t (s)-n oo (s)|. 

The mixing time T m i x (e) is defined by for any e > 0, 

T mix (e) := min{t|A(t') < e for all t' > t}. 

In order to study an upper bound and a lower bound 
of the mixing time, we introduce a relaxation time of 
a Markov chain. Denote the eigenvalues of a transition 
matrix V by \ (i = 1, 2, • • • ) in decreasing order. When a 
transition matrix is irreducible and aperiodic, it is known 
that 1 = Ai > A2. A relaxation time T re i is defined by 

T rc i = (1 - max IAjI) -1 , 

which gives bounds of the mixing time T m ; x (e) such that 

(T re] - 1) log(^) < T mix (e) < log(— —)T Ieh (CI) 
ze eii m i n 

where II m i n := min s IIo^s) is the minimum stationary 
distribution j5fj| . 

Although the relaxation time provides both of the up- 
per and the lower bounds of the mixing time, it does 
not give tight bounds. Hence, we introduce a coupling 
method for investigating the upper bound of the mixing 
time. A pair of two random walks (X t ,Y t ), where t de- 
notes the number of steps, is said to be a coupling of a 
Markov chain when the following two conditions are sat- 
isfied. First, each of X t and Y t is a faithful copy of the 
Markov chain. Second, (X t , Y t ) should satisfy the condi- 
tion that X t = Y t implies X t +\ = Y t +\. For a coupling 
(X t , Y t ), we define a stopping time T xy by 



Tr 



mm{t\X t = Y t , when X Q = x, Y = y}. 
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By definition, X t = Y t for all t > T xy . The stopping 
time is related to the mixing time through the following 
theorem (50l |. 

Theorem 4 Let (X t ,Y t ) be a coupling of a Markov 
chain and T xy be the stopping time. Then, 

A(i) < maxProbp 1 ™ > t]. 

Since the mixing time is obtained from A(t), we can de- 
rive an upper bound of the mixing time from the stopping 
time. 

In the following, we use a relaxation time to obtain a 
lower bound of the mixing time and investigate an upper 
bound of the mixing time by using the coupling method. 

Appendix D: Coupon collecting starting with 
non-zero coupons 

We consider a coupon collecting and show that for c > 

0, 

Prob[r c ^) n > r(bg(r - r ) + c)] < e" c . (Dl) 



Suppose that any of the j-th coupons (j = 1 , ■ • • , r — 
r ) are initially not at hand. We denote by Aj an event 
where the j-th coupon has not been collected within 
r(log(r — ro) + c) steps. An upper bound of the left hand 
side of Eq. (|D1[) is obtained by 



Prob[r c ^ n >r(log(r-r ) + c)] 
= Prob[U^Z[ ^] 

r— rp 

3=1 
r— rp 

3=1 T 

= (r — ro)(l — iy r ( lo g( r - r o)+ c ) 
r 

< (r - r ) exp[- log(r - r ) - c] 
= e" c . 
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